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Abstract 

A wrong derivation of the phase of a propagating massive particle 
which has repeatedly appeared during the last years has the advantage 
of leading once more to the important question of the phase difference 
between the mass eigenstates constitutive of an oscillating neutrino de- 
scribed in the plane wave formalism. Serious errors of principle are pointed 
at in a number of simple calculations and the oscillation length is derived 
in a way which allows to show that the standard result can suffer but 
minute variations in the ultrarelativistic case. 
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1 Introduction 



Neutrino oscillations now seem to be firmly established on experimental grounds. 
It is therefore of necessity to know how to extract physical parameters, viz 
masses and mixing matrix elements, from the data. Of paramount importance 
for this purpose is the formula linking the measured oscillation parameters to 
the unknown masses of the eigenstates constitutive of an oscillating neutrino. 
This famous 'oscillation length formula' has been challenged by many authors 
and its derivation is often bungled by others. The purpose of this paper is 
not to challenge the 'orthodox' point of vue, but to try to derive the oscillation 
phase as neatly as possible in the plane- wave formalism, baring hand- waving and 
faulty approximations. It is not a difficult exercise, but it is, curiously enough, 
at variance with a number of papers by authors who, knowing the result in 
advance, manage to find the right answer for the wrong reasons, doing errors 
which arc often the same as those of the proponents of alternative formulae and 
are based on the same ill employed approximations. Prominent among these are 
the so-called ultra-relativistic approximations which amount, for the careless, 
to the mere replacement of v by c = 1 without any consideration for what they 
were taught about doing asymptotic expansions. The first purpose of this paper 
is therefore pedagogical. It is deeply deplorable that so many wrong derivations 
are still plaguing seminars, articles and books on this question, especially when 
papers appear in journals devoted to the teaching of physics. 
However, this article is also aimed at physics itself. We hold it that doing proper 
maths helps to understand the problem. That the phase is not what the authors 
of wrong papers want it to be forces to employ a little more thought; there 
is always something to be learned from such exercices. With very elementary 
maths, we end up showing that there is no room for deviations from the standard 
oscillation length with light, ultrarelativistic neutrino-like particles, that is, in 
the only case where an oscillation length has a physical meaning. 

2 Phase of a definite mass particle 

The derivation of the phase of a propagating massive neutrino published in [1] 
(see also [2, 3, 4, 5, 6, 7, 8, 9, 10] and the book [11] p.290 eqn 8.37 and probably 
many more which reproduce the same calculation) is plainly wrong. 

Starting from the relative phase 1 

$(£,T) =pL- ET (1) 

acquired by the state vector of a particle of energy E and momentum p upon 
displacement from space-time point (0, 0) to space-time point (T, L), the author 
makes the 'ultra-relativistic approximation' T = Lon the ground that neutrinos 
are ultra-relativistic particles, and from there, using the dispersion relation E 2 = 

1 A c = H = 1 system of units is used 
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p 2 +m 2 andp f=a E for an ultrarelativistic particle, he gets the following equation: 
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*=(p-^ = -— (a) 

We claim that pL — ET = (p — E)L is an ugly equality which should 
be banned from decent physics littcrature although it can be found in many 
places, starting with the review dedicated to neutrino masses and oscillations 
in various editions of the celebrated Review of Particle Physics (see eg [9], [10] ) 

It is quite clear that this cannot be correct: 



• First and foremost, L is vT, not T and replacing the variable by its limit 
in calculating an asymptotic expansion is not allowed. If he were self- 
consistent, the author of [1] would do the same with p = vE, afterwards 
the phase would collapse.. 

• Second, the result must fit seamlessly with what is obtained without the 
ultra-relativistic approximation. But then, when v < 1, anybody would 
use this obvious L = vT which, when combined with the p = vE relation 
also valid for any velocity, yields 

9 v —m 2 T —m 2 L , . 

$ = ET(v 2 - 1) = — — = — — (2) 

This is equal to twice the expression found in [1] (see eqn. (a)) in the 
v = 1 limit, and there is nothing to tell us when, and for which bizarre 
reason, one should suddenly drop this factor when v gets close to one. 

• Third, Lorentz invariance dictates precisely what the phase of a plane wave 
must be in an arbitrary frame. Any wave equation yields the Schrodingcr 
time dependence exp(-imT) in the rest frame 2 and consequently the 
phase in any frame is given by (2), not by half of it. Here L is understood 
as the spatial distance traveled by a classical particle of energy E and 
linear momentum p in a time T. We shall give below a proof of this 
Lorentz invariance which allows to demonstrate (without approximations) 
the point that the author of [1] wanted to make in his paper; but it must 
be stressed that Lorentz invariance of the phase is a logical requirement 
as basic as invariance of the numbers of nuclei in a given lump of matter. 
Although the author of [1] claims that (1) is Lorentz invariant 3 , he does 
not seem to mind the fact that his 'approximation' yields only half the 
Schrodinger phase in the rest frame. 



2 r is the proper time 

3 Of course, Lorentz invariance of (1) has no meaning as long as L and T or E and p are 
not somehow related. 
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Curiously, a similar treatment employed in the Review of Particle Physics 
has been challenged, not because of its physical and mathematical unsoundness, 
but on the ground that the 'overall' phase is not observable [12] and is not 
gauge-invariant and that readers were induced by the text into taking "the 
phase seriously". It is clear however, that this is not an 'overall' phase, but 
the phase difference acquired in a specific space-time translation and that it 
is precisely this phase which is used when building a wave packet or deriving 
a space-time propagator from the path integral. Since it is a difference, it is 
obviously gauge-invariant (whatever the relevance of this concept for a free 
non-interacting neutrino..) and should be "taken seriously" as it is done in, e.g., 
any derivation of the solutions of a wave equation starting from the rest frame. 
Anyway, even non gauge-invariant constructs such as the electromagnetic 
potentials have to be taken seriously as far as basic maths are concerned. Their 
inobservability does not justify hand-waving calculations. In the case at hand, 
we all know that there is no such thing as an intrinsic (rather than " overall" ) 
phase of the state vector, but phase variations do have a meaning and cannot 
be multiplied by arbitrary factors. 

Observe also that no justification is givent in [1] for the use of an 'ul- 
trarelativistic approximation' (v=l) in the relation between elapsed time 
and distance and its not being used in the relation between energy and 
linear momentum. The reason for using the first approximation is evidently 
that it later allows to calculate the phase difference (between two mass 
eigenstates) at a unique same space-time point without further discussion; 
the unstated reason for not using the second approximation is that the 
phase would collapse to as for light rays, if both were used simultaneously. 
(This second approximation is in fact also used, but only in the very last 
step of the derivation, when p is safely added to E instead of being subtracted..) 

Therefore, the reason for the difference and the factor of two boils down to 
1 — v 2 f=a 2(1 — v) when v — ► 1. Using only one factor of v amounts to dropping 
the 2. 

The same trick is used with the same lack of care and comments in many 
such 'simple' derivations of the oscillation length. Physicists should know that 
an equivalence is not an equality, that a statement about a limit is a particular 
kind of equivalence and that even a physical quantity never becomes equal to 
its limit under any approximation when this limit is outside the physical range. 
'Equal in the ultrarclativistic regime' is devoid of meaning and fraught with 
traps. 

3 Lorentz invariance of the phase 

Note, however, that we do not challenge the main point of [1] viz. the 
Lorentz invariance of the oscillation probability, which must hold for obvious 
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consistency reasons. We shall start by showing, along the lines of [1] but 
with a little more care and without any 'ultra-relativistic approximation', 
that the Lorentz invariant quantity is really T/E = L/p where (L,T) are 
the space-time coordinate differences between the events of production and 
detection. Contrary to the claim of [1], his 'approximation' T = L cannot be 
the crucial point in a check of Lorentz invariance, which must also hold between 
two frames wherein the particle velocities are not close to c 

Let us assume that the object is produced at point x = and time t = 
and that its speed is V in frame K where the stationnary target is situated at 
x = L. Let T = L/V be the impact time in K. 

In a frame K' with axes parallel to those of K and moving with velocity v 
along Ox, with the origins so choosen that t' = t = 0, the impact time will 
be V = = T -^0=T 0n thc othcr hand , thc object energy in K' will 

be E' = E ~ vp . i = E if its energy and momentum are E and p in K. 

Therefore T' / E' = T/E = L/p = V /p' 4 as announced and this immediately 
implies Lorentz invariance of the phase as shown by formula (2). The true 
invariant is T/E = L/p and not, as stated in [1], L/E to which it is equivalent 
(not equal !) in the ultrarelativistic regime. For further reference, observe that 
the equality holds because the same V is used in both position and momentum 
space. 

4 What is the "equation of the trajectory" ? 

We are dealing with quantum objects and one might criticize our use of v = p/E 
in the time-distance relation; however, it is well known that, when one uses 
wave packets instead of plane waves, the center of the free wave packet follows 
the classical trajectory, this being true both in non-relativistic and relativistic 
quantum mechanics 5 . The corresponding velocity is thc group velocity of thc 
packet as can be shown by a stationnary phase argument or by calculating the 
expectation value of the position operator. Using v = = ^ in the plane 
wave approximation where the time-distance relation must be added by hand 
is therefore well justified. Hence, if a phase is to be used for the center of thc 
wave packet describing a one particle state, it cannot be anything else than (2). 
But we shall see that it is quite irrelevant for the oscillation problem. 

4 The elementary derivation just given can be mimicked for L and p. The second equality 
trivially follows from p = vE and L = vT 

5 In the latter case, and depending on how the wave- packet is built, there is a quite negligible 
additional oscillation around the central value, (zittcrbewegung), with an amplitude of thc 
order of h/E that is, 2fm for a 10 MeV particle 
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5 Phase difference of different mass eigenstates 



5.1 Standard approximations 

Standard derivations of the oscillation formula in the plane-wave formalism re- 
sort most often to approximations known as "equal E" or "equal p" in which 
one postulates equal energies or equal momenta for the two mass-cigcnstatcs 
constitutive of the flavour neutrino, notwistanding the contradiction with pro- 
duction kinematics and the fact that these hypotheses are quite arbitrary - 
these equalities aren't even Lorentz invariant. The reason for so doing is that 
a connection between time and space must be postulated since it cannot be 
inferred in this case from a dynamical calculation. With the above mentionncd 
approximations, the spatial or the temporal interval disappears upon subtract- 
ing the phases of the two mass eigenstates and there is no discussion about 
the (space-time) point where one chooses to calculate the interference pattern. 
Specification of the distance is enough. Using his above equation (a) in the 
so-called " equal E approximation" , the author of [1] immediately falls back on 
the standard result that gives rise to the neutrino oscillation length quoted in 
hundreds of papers on the subject. There are, however, less qucstionnable ways 
of retrieving the same result. 

5.2 What is wrong with these phases ? 

As we have shown, the phase is given by (2) and at first sight, subtracting 
similar expressions written for m = mi and m = m-x will yield twice the 
standard result for the oscillation phase. 

To get their factor straight, the authors of [13] introduce a "space velocity v s " 
different from the phase velocity and write <f> = — ET(l/j 2 + l/j^)/2 after 
what 7 S , which is the same for all mass eigenstates, cancels in the substraction, 
yielding the requested factor 1/2. Indeed, v s = L/T is the same for all mass 
eigenstates "as required by quantum mechanics" write the authors of [13]. 

We prefer a less contorted approach, using directly the already established 
results. A same "space velocity" would preclude the separation of the centers 
of the wave packets and the washing out of oscillations which will necessarily 
occur in a complete treatment of the problem. 

However, sticking to what we have so far established, the result (2) does not 
mean that the standard oscillation phase is wrong and that proponents of the 
formula modified by the famous factor of 2 are right. The point is that for 
different values of m (and a given observation point, that is, a given I in a 
frame were the target is stationnary) the factor T in (2) cannot be the same. 
Use has been made of the kinematics (T = L/v) of the center of the wave packet 
and since the energies are supposedly equal, different masses lead to different 
velocities, hence to different travel times. Likewise, using the second form of 
(2), the equal L factors hide different times because of different velocities. 
But phases must be compared at the same space-time point. It is therefore 
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meaningless to subtract them after having reduced them to the form (2). 

It is quite obvious that there is no way to derive any interference from 
classical kincmatical calculations. Classical particles are concentrated, so to say, 
at the center of their wave packets and therefore cannot interfere as soon as these 
centers separate. On the other hand, what is required by quantum mechanics 
is to calculate overlap amplitudes at a given point and a given time. Using 
different instants and/or different space points simply makes no sense. This 
is the true justification for the same "space velocity" of [13] and this is where 
most proponents of alternative formulae stumble. They forget that elementary 
particles have little to do with sand grains and that it is only the centers of 
the wave packets which follow L = vT trajectories. But it is precisely that 
" fuzincss" which allows to calculate a phase difference and to find interferences 
and oscillations. Classical particles cannot interfere. 



5.3 Various solutions 

The solution here is to go back to the original expression (1) of the phase where 
L and T can be fixed independently. Then, in the equal E approximation, the 
time-energy parts cancel and the space-momentum parts yield the answer: 

S^ = L P1 -L P2 - L{P "- P ^ rmi-mj 



Pi + Pi 2p 

where use has been made of the 'equal E' hypothesis and an obviously defined 
p has been introduced. Note that we have not used any relativistic approxi- 
mation to expand the momenta. The expression found is valid in principle for 
hypothetical non relativistic neutrinos ([14]) 6 . 

The same reasoning applies to the "same p" approximation. Formula (2) 
written in terms of p and L does not lend itself to a simple subtraction between 
mass eigenstates, because the two subtracted expressions do not correspond 
to the same time for the already explained reasons. Here again, one must go 
back to the original expression of the phase, whereoff the pL parts cancel upon 
subtraction to yield: 

El - E? mi - m? 
5$ = TE 2 - TE\ = T-2 -2- = T- 



Ex + E 2 2E 

The last expression must be converted to a function of distance to make connec- 
tion with experiments. Using the average velocity defined by p/E, one retrieves 
the expression found above with p — > p, again without a relativistic expansion, 
(see [14]) 

Even the (notoriously unphysical, see [15]) "same velocity" approximation 
yields the standard answer when treated with care [14]. The reason for which 
a wrong hypothesis leads to a reasonable result will become clearer in the next 
section. 



3 See however the next section for a caveat 
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6 A better Ansatz 



In [14, 16] the present author (unaware of the existence of [17]) showed that not 
only could one retrieve the standard result without any relativistic expansion 
as shown here above, but also without any of the "equal energy" or "equal 
momentum" or "equal velocity" approximations, by postulating an appropriate 
velocity for the center of would-be wave packet. This velocity is defined as 
v ° = e\ ~+e-> ( see a ^ so anc ^ [1^] wnere this same value is used to disprove the 
existence of oscillations for the recoil particle) It was guessed on the basis of the 
average velocity found as already described in the equal p approximation, or 
in the equal v approximation were it is derived from the (by itself unrealistic) 
assumption that |r- = |J 7 

As already said, the usual equality assumptions contradict production kine- 
matics, be it through a reaction or through a decay. No such violence is done 
to kinematics in the present case. Energy and momenta can be given their cor- 
rect, conservation laws respecting values. The particular expression for vq arises 
quite naturally from writing the phase difference as: 

Sp 2 SE 2 
<5$ = xSp - tSE = x-^r - t—— 
Yip IjE 

Converting this to a function of x alone, one is led to write i = § and the choice 
v = |^ tends to impose itself, yielding: 

dm 2 5m 2 
S * = ~ X ^p~ = - X ^f 

which is the standard oscillation phase with the bonus of showing that 
the kinematical variable to use, would the distinction between energy and 
momentum be relevant, should be the latter, since this calculation is valid for 
any velocities ([20], [14, 16]) 8 . 

Although natural in the above context, the value Vq springs out of the blue 
and is not justified by any convincing calculation. It should come out of the 
dynamical evolution equation in a wave packet or field-theoretical treatment, 
but since our purpose is to keep things simple, we will not embark on such 
calculations and refer the reader to the enormous existing litterature which we 
don't even feel able to review properly. For a very complete report of elaborate 



7 The reason for which this hypothesis is wrong in principle is that strictly equal velocities 
entail masses in the same ratio as energies that is, almost equal for a beam of more or 
less defined energy - and equal masses are obviously completely unrealistic (see [15]). The 
practical reason for which it works is that the inverse function 7— > v = \f\ — I/7 2 is fiat and 
undistinguishably close to 1 well before 7 reaches values which are relevant to the neutrino 
case. Therefore very different E/m's yield the 'same' (for all practical purposes) velocities 

8 Note however, that an oscillation probability independent of the production and detection 
processes is an illusion for hypothetical neutrino-like particles of sizeable mass. No factoriza- 
tion of production, propagation and detection probabilities can hold here 
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treatments, see [21]. Here, we will content ourselves with showing that devi- 
ations from this 'desirable' velocity, as Lev Okun calls it, do not lead to any 
serious alterations of the standard oscillation length . 



7 Oscillation length 

vo as above defined can trivialy be looked at as an average of the classical 
velocities of the two mass eigenstates weighted by their energies: 

Ei E 2 

Vq = V\ H V2 

Ei + E2 Ei + Ei 

with Vi = j^j-. Hence it lies between these two, like any reasonable value for the 
classical velocity of a mixed state. We shall therefore parametrize an arbitrary 
(but 'reasonable') velocity by 1/v = l/v + e which yields the phase shift as a 
function of distance : 

5m 2 SE 2 Sm 2 SE 2 

d<±> = —x x e = — x x wn£ 

2p EE 2p 2p 

yielding an oscillation length: 

Apm Apn 



L = 



5m 2 + vo5E 2 e 5m 2 + 2p8Ee 
which differs from the standard result by the second term in the denominator. 



Assume now that the neutrino is produced in a reaction with total c.o.m. 
energy squared s. Then Ei = (s + mf — m 2 )/{2y/s) where m r is the effective 
mass of the system recoiling against the neutrino. Then SE = Sm 2 j '(2-y/s) and 
we have: 

5m 2 (l + pe/ x /s) 1+pe/x/s 
where Lo is the standard oscillation length . 

Now for light neutrinos, p w E « l/2(y / s — m 2 /^/s); p/y/s is less than 1/2 
and can be fairly lower. For example, in a ir — * \xv decay, it is about .21 and it 
can be still much lower in three body decays. 

|e| itself can be easily bounded from its definition. But it clearly can't exceed 
\v\ — V2\j{viV'2) which is extremely low since with standard neutrino masses, 
many figures will cancel between the two w's in the numerator. Therefore there 
isn't much room. Only very small numerical variations w.r.t. the standard 
oscillation length value are possible in this frame 9 . One sees that the reason 

9 At 10 GeV c.o.m. energy, a v assumed to have an 19 MeV mass (limit derived for 
i/ T 's from r decays but enormously larger than expected standard neutrino masses given the 
measured Sm? 's and the limits on m Ue ) recoiling against a proton has v = ,99999a; where x 
stands for the first digit not equal to 9 



which makes this oscillation length so sharply defined is the same that renders 
the equal velocity hypothesis practically useable in spite of its being wrong in 
principle. Indeed, v = 1 after all, even if it is not permitted to write it before 
the end of the calculation ! 

8 Lorentz invariance of the oscillation phase 

Be it calculated with the usual equality hypotheses or with the ansatz proposed 
in section 6 for the classical velocity of the mixed neutrino, the oscillation phase 

c 2 

is: (5$ = —x^- and it is Lorentz invariant, as shown in section 3, as soon as x 
is interpreted as the distance traveled by the neutrino between production and 
detection in the frame where it has momentum p. It has to be so since this phase 
directly gives the probability that an interacting neutrino yields a lepton of the 
initial flavour. The 'oscillation length' L OSCl derived by equating this phase to 
2tt has the meaning of the source-detector distance giving a maximum of this 
probability only in the frame where the target is stationnary. In other frames it 
can be found by scaling L osc in the same ratio as the momenta and it represents 
the distance traveled by the neutrino during its existence. 

9 Conclusion 

Having recently assisted to a seminar where the untenable derivation of the 
oscillation phase as above described was proposed, the present author has 
decided to write this note in the hope of contributing to bar widespread and 
very unpedagogical hand-waving in the simple way of attacking the question of 
vacuum neutrino oscillations. He took the occasion to show why only minute 
deviations from the standard oscillation phase are possible in this simplified 
frame and to straighten up the derivation of the Lorentz-invariance of the 
oscillation probability. 
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